Formation and rotation of skyrmion crystal in the chiral-lattice insulator Cii20Se03 
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(Dated:) 

Small angle neutron scattering experiments were performed on a bulk single crystal of chiral- 
lattice multiferroic insulator Cu20Se03. In the absence of an external magnetic field, helical spin 
order with magnetic modulation vector q || (001) was identified. When a magnetic field is applied, 
a triple-^ magnetic structure emerges normal to the field in the A-phase just below the magnetic 
ordering temperature T c , which suggests the formation of a triangular lattice of skyrmions. Notably, 
the favorable g-direction in the A-phase changes from q || (110) to q || (001) upon approaching T c . 
Near the phase boundary between these two states, the external magnetic field induces a 30° -rotation 
of the skyrmion lattice. This suggests a delicate balance between the magnetic anisotropy and the 
spin texture near T c , such that even a small perturbation significantly affects the ordering pattern 
of the skyrmions. 



PACS numbers: 75.25.-j, 75.85.+t, 74.25.Ha 

The skyrmion is a topologically stable particle-like ob- 
ject, which appears as a vortex-like swirling spin texture 
in chiral-lattice magnets. 1 — In metallic compounds, con- 
duction electrons flowing over skyrmion spin texture feel 
a fictitious magnetic flux arising from the spin chirality of 
the skyrmion or skyrmion number, thereby producing the 
topological Hall effect. 4,5 On the other hand, the electric 
current can drag the skyrmion via a spin transfer torque, 
while the motion of the skyrmion generates a transverse 




FIG. 1: (Color Online) (a) Crystal structure of Cu 2 OSe0 3 
with two distinct Cu 2+ sites with different oxygen coordina- 
tion, (b) Proper screw spin order with a single magnetic mod- 
ulation vector q. (c) Skyrmion crystal with triple magnetic 
modulation vectors <?i , <?2 , and qs , which lies within a plane 
normal to the applied magnetic field (H). The background 
color represents the out-of-plane component of the local mag- 
netic moment (m z ). 



electromotive force as the emergent electric fieldAi Such 
electric controllability, as well as its particle-like nature 
with nano metric size, highlights skyrmion as a unique 
magnetic object with potential applications in spintron- 
ics and high-density magnetic storage devices. In chiral- 
lattice magnets, the magnetic interaction acting on two 
neighboring spins (Si and Sj) is mainly composed of two 
terms; the Si • Sj-like (ferromagnetic) exchange interac- 
tion and the Si x Sj-like Dzyaloshinskii-Moriya interac- 
tion. The latter term tends to orient the neighboring 
spins in a canted fashion, and hence the helical spin struc- 
ture with single modulation vector q is generally favored 
(Fig. 1(b)). In a limited temperature (T) and mag- 
netic field (H) window, however, skyrmion spin state can 
dominate over the helical ground state. In most cases, 
skyrmions crystallize to form a triangular lattice within 
the plane normal to the magnetic field (Fig. 1(c)), while 
they can also exist as independent particles^ 

Magnetic skyrmions have been first experimentally 
discovered in intermetallic compounds such as MnSi, 8 
Fei-^Co^Si^^ and FeGeJ^ These compounds share the 
same B20 structure with chiral cubic symmetry, and 
commonly show the helimagnetic gound state. In the 
bulk form, their magnetic phase diagram includes an H- 
induced narrow region, called "A-phase", just below the 
magnetic ordering temperature T c . 12 Small angle neutron 
scattering (SANS) studies in the A-phase revealed that 
the three magnetic modulation vectors form an angle of 
120° with respect to each other within a plane normal to 
the applied H, which suggests the formation of a trian- 
gular lattice of skyrmions (Fig. 1(c)) For thin films, 
such a skyrmion crystal (SkX) has indeed been observed 
in real-space through Lorentz transmission electron mi- 
croscopy (TEM), while the SkX in the thin- film form is 
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FIG. 2: (Color Online) Typical SANS profiles for the (110) 
plane of bulk Cu20Se03, measured at various magnetic fields 
and temperatures. All the data were obtained in the H- 
increasing scans, (a)-(d), (e)-(g), (h)-(i), and (j) represent 
the ones at 56.75 K, 57.25 K, 58.00 K, and 57.75 K, respec- 
tively. The inset shows the measurement setup; fa n and k ou t 
are the incident and scattered neutron wave vectors, respec- 
tively. The magnetic field (H) is applied along fa n || (HO), 
i.e. perpendicular to the observed plane. 



stabilized over much wider T- and i^-range than the bulk 

Q 1 1 

case**i±± 

Recently, through Lorentz TEM imaging of a thin 
film, a SkX has also been observed for the insulat- 
ing compound Cu2 0Se03. 13 Its crystal structure be- 
longs to the same chiral cubic space group P2\3 as the 
B20 compounds, but with a very different atomic ar- 
rangement (Fig. 1(a)) j 14 i 15 The magnetic phase dia- 
gram of Cu2 0Se03, like that of the B20 compounds, 
is radically different between bulk and thin film. For 
bulk specimens, previous powder neutron scattering and 
NMR studies suggested a three-up one-down type fer- 
rimagnetic arrangement between two inequivalent Cu 2+ 
(S = 1/2) sites^"— while ac magnetic susceptibility 
measurements indicated a more complex magnetic phase 
diagram characterized by an A-phase similar to that of 
the B20 compounds^ 3 - The magnetic structure in the 
A-phase of bulk Cu20Se03 has not yet been resolved, 
although the analogy with the B20 compounds implies 
the possible formation of a SkX in this region. Notably, 
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FIG. 3: (Color Online) Magnetic field (H || [110]) dependence 
of (a) ac magnetic susceptibility x an d (b) integrated inten- 
sity / of magnetic reflections obtained from SANS profiles. 
Temperature dependence of x an d / is also plotted in (c) 
and (d). 



the bulk form of Cu20Se03 exhibits a magnetically in- 
duced electric polarization in the A-phase ^ 13 ! 16 and thus 
is considered as a multiferroic (i.e. material with both 
magnetic and dielectric orders) J^r— Such magnetically- 
induced ferroelectrics often allow electric field control of 
magnetism^"— and the observed magnetoelectric cou- 
pling in the A-phase implies the possible manipulation 
of skyrmions by external electric fields in insulators. To 
establish the coupling between electric polarization and 
skyrmions, a detailed investigation of spin textures in 
bulk Cu20Se03 is urgently required. 

In this paper, we study the long- wave length magnetic 
structure of bulk Cu20Se03 by SANS measurements. 
We identified the helical nature of the magnetic ground 
state, and further discovered the formation of a skyrmion 
lattice in the bulk A-phase. This establishes the coupling 
between skyrmion spin texture and electric polarization, 
which promises the potential manipulation of skyrmions 
by external electric fields, not electric currents. Interest- 
ingly, the skyrmion lattice shows a 30°-rotation within 
the plane of observation as a function of temperature and 
magnetic field. This finding suggests an extreme sensi- 
tivity of the skyrmion ordering pattern to perturbation. 

Single crystals of Cu20Se03 were grown by the chem- 
ical vapor transport method^ AC magnetic susceptibil- 
ity x' ( a t 700 Hz) was measured with a SQUID mag- 
netometer. SANS measurements were performed at the 
diffractometer MIRA at FRM II using neutrons with a 
wavelength A ~ 5.1 A. Backgrounds were determined at 
high temperatures and subtracted accordingly. A mag- 
netic field was applied parallel to incident neutron beam 
using Helmholtz coils (for setup see Ref. 8). Data were 
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FIG. 4: (Color Online) Temperature (T) and magnetic field 
(H) phase diagram obtained for the bulk Cu20Se03 with 
H || [110]. Diamonds represent the anomalies in ^-profile, 
and circles indicate the data points where SANS measure- 
ments were performed. In both cases, closed and open sym- 
bols are obtained from the H- and T-scans, respectively. The 
background color indicates the properties of SANS profile at 
each data point; color saturation represents the integrated in- 
tensity of six- fold magnetic reflections, and the hue stands for 
the rotation angle (0) of the magnetic Bragg position with 
respect to the (110) direction (see the bottom panel). Typi- 
cal SANS profiles at selected data points are shown in Figs. 
2 (a)-(j). In the data points with cross symbols, the mixed 
SANS pattern as shown in Fig. 2 (f) was observed. Here, 
the A-phase consists of SkX(l) and SkX(2); The former cor- 
responds to the skyrmion crystal state favoring q || (110), and 
the latter does the one favoring q || (001). 



recorded for fixed sample orientations, following rocking 
scans with respect to a vertical axis of typically ±1.8°. 
All measurements were performed on the same as-grown 
crystal with the size of 8 mm x 5 mm x 3 mm, and the 
magnetic field was always applied normal to the (110) 
widest surface. All i^-scans were performed in the H- 
increasing process after zero field cooling (ZFC), and the 
T-scans were performed in the warming process after 
ZFC and subsequent application of H. 

First, we focus on the behavior at 56.75 K, which is 



close to T c ~ 58 K. Figure 2(d) indicates the SANS profile 
in zero magnetic field. We found two magnetic Bragg re- 
flections with q || (001), which suggests that the magnetic 
ground state of bulk Cu2 0Se03 is not simply ferrimag- 
netic, but rather helically modulated. The observed wave 
number q ~ 0.010 A -1 corresponds to the modulation pe- 
riod of ~ 63 nm, which roughly agrees with the helical 
modulation period ~ 50 nm previously reported for the 
thin film J* 3 ** 7 This value is much longer than the Cu-Cu 
atomic distance and thus compatible with a local ferri- 
magnetic arrangement. The direction and wave number 
of q remain unchanged down to 10 K (not shown). Due to 
the high symmetry of the cubic lattice, the equal popula- 
tion of q || (001) domains will coexist. Notably, the prop- 
agation vector of the helimagnetic ground state is differ- 
ent between bulk (q || (001)) and thin film (q || (110)) 13 
for Cu2 0Se03. Recent theories suggest that the stable 
g-direction in the helical ground state depends on the 
dimensionality of the system^ They predict q || (110) 
for a two-dimensional (2D) system, while q || (001) or 
q || (111) for a 3D system. The observed q directions 
in the bulk (3D) and thin film (2D) forms of Cu 2 OSe0 3 
agree well with these predictions. 

In the following, we discuss the development of SANS 
profiles under magnetic field applied along the out-of- 
plane (110) direction. Previous Lorentz TEM mea- 
surements on a thin film suggested a screw-type helical 
ground state, where the spins rotate within a plane per- 
pendicular to q (Fig. 1(b)). Since antiferromagnetically 
aligned spins tend to orient normal to H, the applied 
magnetic field should favor the H \\ q relationship and 
reorient q along the out-of-plane direction in the present 
SANS setup. Figure 3 (b) indicates the independence 
of the integrated intensity of the magnetic reflections ob- 
tained from the SANS profiles. Above 50 Oe, the mag- 
netic reflection vanishes from the plane of observation 
as shown in Fig. 2(c), which is consistent with such an 
^-induced reorientation of q in the helical state. 

Upon increasing the magnetic field further, magnetic 
reflections in the SANS profile emerge again at 200 Oe 
< H < 300 Oe (Fig. 3(b)). The corresponding in- 
dependence of the ac magnetic susceptibility x' shows a 
clear dip anomaly in this field region (Fig. 3 (a)), which 
signals the transition into the A-phase j 13 i 3Q Figures 2 (a) 
and (b) indicate the SANS profiles at 255 Oe and 213 Oe, 
respectively. We observed a six-fold symmetric pattern 
of magnetic Bragg spots within a plane normal to the 
applied which evidences the formation of a skyrmion 
crystal in the A-phase of bulk Cu 2 OSe0 3 (Fig. 1 (c))&^ 
Three magnetic modulation vectors form an angle of 120° 
with respect to each other, and one of them satisfies the 
q || (110) relationship. Interestingly, the favorable q di- 
rection in the SkX state of bulk Cu20Se03 depends on 
temperature. Figure 2 (j) shows the SANS profile at 
T = 57.75 K with H = 170 Oe. The six-fold magnetic 
reflections verify the skyrmion lattice formation, while 
the SANS pattern rotates by 30° with respect to the one 
at 56.75 K. In this case, one of the magnetic modulation 
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vectors rather favors the q || (001) relationship. 

Based on the various T- and iTscans of x' an d the 
SANS profiles, we have summarized the H-T phase dia- 
gram for bulk Cu 2 OSe0 3 with H || (110) (Fig. 4). The 
anomalies in the x' profile are indicated with diamonds, 
and the A-phase is identified in the narrow region of 56 
K < T < 58 K and 100 Oe < H < 350 Oe. The back- 
ground color in the phase diagram indicates the nature 
of the SANS profiles at each data point. The color sat- 
uration represents the integrated intensity of the six-fold 
magnetic reflections, and the hue reflects the rotation 
angle of the magnetic modulation vector with respect to 
the (110) direction (see the bottom panel of Fig. 4). The 
emergence of six-fold magnetic reflections always coin- 
cides with the dip anomaly in the x'-profiles (Figs. 3 
(a)-(d)), which identifies the A-phase as the SkX state in 
the bulk crystal. Our SANS measurements have further 
revealed that the A-phase is divided into two regions; the 
SkX(l) phase favoring q || (110), and the SkX(2) phase 
favoring q || (001). Near the phase boundary, we can 
even switch these two SkX states by the external mag- 
netic field. Figures 2 (e)-(g) indicate the independence 
of the SANS profiles at T = 57.25 K. The transition from 
SkX(2) to SkX(l) is induced as H increases from 170 Oe 
to 255 Oe, and a mixed pattern of these two profiles due 
to phase coexistence is observed at the phase boundary 
(213 Oe). Note that neither x' nor the integrated in- 
tensity of the magnetic reflections show any noticeable 
anomaly at the transition between these two SkX states 
(Figs. 3(a) and (b)). The magnetic modulation period 
remains unchanged among the helical and the two SkX 
spin states. 

For metallic MnSi, a rotation of the skyrmion lattice 
by up to 10° has recently been reported under the com- 
bination of a temperature gradient AT and an electric 
current j. 6 Its rotation angle becomes larger as the mag- 
nitude of AT and j increases, and the origin of this phe- 
nomenon is ascribed to a Berry phase collection by con- 
duction electrons flowing over skyrmions. This scheme is 
clearly inapplicable to the present case, since Cu20Se03 
is an insulator and the rotation angle seems to be fixed 
to 30° without any current. According to a recent phe- 
nomenological Ginzburg-Landau analysis, the pinning of 
the (/-direction in the SkX state is governed by a mag- 
netic anisotropy of sixth order in spin-orbit coupling£ii£ 
For example, the terms like Y^qi^x + Qy + Qt)\™q\ 2 an d 

E g (qii 2 y + <iUl + ai<il)\™ q \ 2 favor q || (110) or q || (001) 
depending on the sign of their prefactor. Here, m q is 
the Fourier component of the spatially dependent local 
magnetization M(r). The observed rotation of SkX sug- 
gests that the pinning of q is very weak, and that the 
SkX(l) and SkX(2) states are nearly degenerate. Their 
relative stability is reversed near T c , perhaps because the 
enhancement of thermal fluctuation 8 and/or higher har- 
monics in spin texture^^i near T c provide nonmonotonic 



contributions to the anisotropy terms. While a further 
theoretical analysis is necessary to fully clarify its ori- 
gin, the current observation of a SkX rotation indicates 
an extreme sensitivity of the skyrmion ordering pattern 
against external field or perturbation. 

Combined with the previous reports for the thin film, 13 
the current discovery of long- wavelength spin textures in 
a bulk crystal fully establishes Cu20Se03 as an insulat- 
ing counterpart of the B20 compounds. Both systems 
have a chiral cubic lattice, and their magnetic phase di- 
agrams are characterized by the helical ground state as 
well as the iTinduced SkX state. The stability of the SkX 
state depends on the sample dimension and it shrinks 
into the narrow A-phase in the bulk limit j 11 1 13 These 
similarities are rather surprising considering the signif- 
icant differences between the two systems, such as the 
electronic structure (metal or insulator), the local spin 
arrangement (ferrimagnetic or ferromagnetic), and the 
atomic arrangement in the crystal. This finding strongly 
suggests that chiral cubic ferro/ferrimagnets may ubiq- 
uitously host skyrmions, irrespective of other details. In 
addition, the identification of skyrmions in the A-phase 
of bulk Cu2 0Se03 establishes the coupling between the 
skyrmion spin texture and the previously reported elec- 
tric polarization, 1 - which promises the manipulation of 
skyrmions by external electric fields. Since an electric 
field in insulators causes negligible joule heat loss com- 
pared to the current-driven approach in metals, the cur- 
rently established magnetoelectric skyrmion may be use- 
ful for the design of unique spintronic devices with high 
energy efficiency. 

In summary, we have performed a detailed small-angle 
neutron scattering study of the bulk form of the chiral- 
lattice multiferroic insulator Cu2 0Se03. The helical na- 
ture of the magnetic ground state, as well as the forma- 
tion of a skyrmion lattice in the A-phase, have been suc- 
cessfully identified. This establishes a coupling between 
electric polarization and skyrmions in the A-phase, which 
implies the possibility of manipulating skyrmions by ex- 
ternal electric fields. Notably, the favorable (^-direction 
in the skyrmion state changes as T c is approached. Near 
the phase boundary, the external magnetic field induces 
a 30°-rotation of the skyrmion crystal. This finding 
demonstrates the delicate balance between the magnetic 
anisotropy and spin texture near T c , and suggests that 
the skyrmion ordering pattern may be controlled by var- 
ious external stimuli. 
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